Correlation between Gamma-Ray bursts and Gravitational Waves 
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The cosmological origin of 7-ray bursts (GRBs) is now commonly accepted and, according to several 
models for the central engine, GRB sources should also emit at the same time gravitational waves 
bursts (GWBs). We have performed two correlation searches between the data of the resonant 
gravitational wave detector AURIGA and GRB arrival times collected in the BATSE 4B catalog. 
No correlation was found and an upper limit /irms < 1-5 X 10 -18 on the averaged amplitude of 
gravitational waves associated with 7-ray bursts has been set for the first time. 

PACS numbers: 04.80.Nn, 98.70.Rz 



I. INTRODUCTION 

Thirty years after their discovery, Gamma Ray Bursts 
(GRBs) are still the most mysterious objects in the Uni- 
verse, as their properties have not yet been associated to 
any well-known object, while for instance pulsars have 
been almost immediately identified as final objects of 
the evolution of stars and quasars have been well framed 
among galaxy nuclei in particular in the Seyfert class. 
Unfortunately "standard" astrophysical objects do not 
exist for GRBs. The "non standard" compact astro- 
physical objects which might reproduce the experimental 
characteristic of GRBs involve black-holes (BH) , neutron 
stars (NS) and massive stars. 

The wide range of characteristics prevents the system- 
atic classification of GRBs. Their energy spectrum is 
continuous and non-thermal, and covers the range be- 
tween 1 keV and 1 MeV; their emission lasts from 10 ms 
to 10 3 s. The Burst And Transient Source Experiment 
(BATSE) has detected GRBs with an event rate of about 
one per day between 1991 and 2000. 

One of the most important results of observations is 
that many GRBs are at cosmological distances |IJ. This 
fact, recently confirmed also by the satellite BeppoSAX 
B, which has detected some optical counterpart with 
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redshifts close to z = 1 [|3|-|6|, implies the emission of a 
huge amount of electromagnetic radiation: in few seconds 
the GRBs release an energy of 10 51 — 10 54 erg. On the 
other hand, this implies that the GRBs are rare events 
in a galaxy: in fact, as the rate of GRB events is of the 
order 1/day, their rate amounts to 1 event/10 6 years per 
galaxy §. 

The theoretical framework which gives a coherent ex- 
planation of the experimental data set of GRBs is the 
so-called "fireball" model: an "inner engine" produces a 
flux of relativistic energy of the order 10 52 erg and it is 
extremely compact and clearly not visible; this relativis- 
tic flux of particles is a kind of "fireball" able to produce 
electromagnetic radiation in a optically thin shell of mat- 
ter. The inner engine works for 1-3 days, producing the 
afterglow, but the bulk of the explosion lasts typically 10 
s (with peaks around 0.5 and 30 s and a variation of six 
orders of magnitude, 10~ 3 4- 10 3 s) when the major part 
of the energy is emitted. 

The inner engine progenitors are likely to be coalesc- 
ing binary systems, such as NS-BH or NS-NS systems, 
or single stars that collapse into BH in supernova-like 
events (the so-called " collapsar" models Q). If this is 
the correct explanation of the GRBs, then we would have 
compact astrophysical systems for which gravity plays an 
important role, and gravitational waves would be emitted 

This scenario has motivated us to investigate the be- 
havior of the gravitational wave detector AURIGA dur- 
ing time spans which include the arrival time of a GRB 
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burst. The AURIGA detector Q is an A15056 reso- 
nant bar of about 2.3 tons with a typical noise temper- 
ature of 7 mK and a bandwidth of about 1 Hz. The 
detector is sensitive to gravitational waves (g.w.) signals 
over 1 Hz bandwidth around each one of its two reso- 
nant frequencies i.e. 913 and 931 Hz. The sensitivity 
of a g.w. detector can be conveniently expressed by the 
quantity h m i ni representing the minimum g.w. ampli- 
tude detectable at SNR = 1, which for the AURIGA 
detector is h min ~ 2 5 x 10" 19 . The AURIGA sen- 
sitivity is enough to detect NS-NS, BH-NS and BH-BH 
mergers which take place within the Galaxy. 

Previous work on GRB-GWB association appeared in 
the literature, specifically about a single GRB trigger 
jy! , and about a set of GRB triggers |l4|-|l6| ; here we 
present an experimental upper limit on such an associa- 
tion obtained with the AURIGA-BATSE data. 

The plan of the paper is as follows: in Sec. II we discuss 
the two methods, coincidence and statistical searches, 
which we use for the association of the GRBs in the 
BATSE catalog with GWBs. The results obtained with 
the two methods are presented in Sec. III. In Sec. IV we 
analize our results and discuss the possibilities opened by 
our methods for future searches. 



II. SEARCH METHODS 

The aim of our analysis is the search of an associa- 
tion between GRB and GWB arrival times within a time 
window W. The analysis has been carried out by means 
of two different procedures which have been already dis- 
cussed in the literature: i) the correlation method JTij] 
which is based on the coincidence between the arrival 
time of candidate gravitational events selected over a 
given threshold and the GRB trigger and ii) a statisti- 
cal method jl6| which relies on a hypothesis test on a 
statistical variable representing the mean energy of the 
gravitational detector at the GRB trigger. 

The coincidence window plays an important role in 
both the analysis. In fact the coincidence window should 
be wide enough to hold the delay distribution between 
GRB and GWB. We notice that its "optimal" value de- 
pends on the astrophysical sources and on the g.w. de- 
tector properties. If we assume that the GRBs are gener- 
ated by internal shocks in the fireball, the delay between 
GRB and GWB is less than 1 second but there are 
still some uncertainties. Moreover, the delay is widen by 
the cosmological redshift. To be as much conservative 
as possible on the distribution of the delays we choose 
W = 5 s. This value turns out to be also consistent 
with the requirements of the filtering procedures of the 
AURIGA detector. In fact, the search for g.w. bursts 
requires a filtering of the data by a Wiener-Kolmogorov 
filter matched to 5-like signals [[l8| and its characteristic 
time is the inverse of the detector bandwidth i.e. ~ 1 s. 
This time, as we shall see below, establishes the timescale 



of the noise correlation of the filtered data. 

In what follows, "<5-like signal" means any g.w. signal 
which shows a nearly flat Fourier transform at the reso- 
nance frequencies of the detector (913 and 931 Hz) over 
a ~ 1 Hz bandwidth. Therefore the metric perturbations 
h(t) sensed by the AURIGA detector can be a large class 
of short signals (of millisecond duration) including the 
latest stable orbits of inspiralling NS-NS or NS-BH, the 
subsequent merging and the final ringdown |l9| ] and the 
collapsar bursts 1 1 which could be expected signals as- 
sociated with GRBs. 



A. Coincidence search 

The coincidence method has been successfully applied 
to the search of coincident excitations of different g.w. 
detectors |2C],^l| ; much work has been devoted to exploit 
its potentialities and to develop robust estimates of the 
background of accidental, even in the presence of non 
stationary event rates jl7|,^l). A g.w. candidate event 
is a local maximum of the filtered data corresponding 
to any excitation of the detector. From the AURIGA 
filtered output we extract a list of candidate events set- 
ting an adaptive threshold in their signal-to-noise ratio 
(SNR = 5). The event lists used in this analysis be- 
long to periods of satisfactory performance of AURIGA 
as described in Ref. 0]. It is worth noticing that the 
AURIGA event search checks each event against the ex- 
pected signal template by means of a \ 2 test (23). The 
event lists contain the information needed to describe a 
5-like signal namely, its time of arrival (in UTC units), 
the amplitude of the Fourier transform, and the detector 
noise level at that time. The BATSE data are taken from 
the 4B catalog by Meegan et al. available on Internet |2~i| ] 
and includes trigger time (in UTC units), right ascension 
and declination, error box and other information about 
every GRB triggered. 

We label with tty the estimate arrival time of the i-th 

candidate g.w. event detected by AURIGA and with 
the trigger time of the fc-th GRB detected by BATSE. 
A coincidence between the i-th AURIGA event and the 
k-th GRB trigger time is observed if — tj\ < W, 
where W is the coincidence window. 

The coincidences found have to be compared with the 
accidental coincidence background due to chance. Two 
standard methods to evaluate the probability of acciden- 
tals have been applied to the pair BATSE-AURIGA: i) 
performing thousands time-shifts of the arrival time of 
one detector with respect to the other and looking for 
accidental coincidences at each shift ]20[ ]; it) assuming 
independent Poisson distribution of event times and us- 
ing the mean measured rates to estimate the accidental 
rate n a 



n a = N A N r 



AT 
T~ 
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where Na is the number of AURIGA events in the period 
T, N 7 the number of GRB events in the same period and 
AT = 2W is the total amplitude of the coincidence win- 
dow (see Appendix |X] for a proof of this basic relation) . 
To avoid the problem of multiple coincidence within the 
same window we prefer to estimate the background with 
W = 1 s and to scale our results with the help of Eq. 
(1) to W = 5 s. The consistency of the two estimates 
ensures that the arrival times of the AURIGA-BATSE 
pair are distributed as Poisson random points even if the 
AURIGA event rate has been found to be not stationary 

0-13 ■ 



B. Statistical search 

The method of the statistical search has been proposed 
for a pair of interferometric detectors by L. S. Finn and 
coworkers |]l6|| ; here we slightly modify their approach to 
the case of a single resonant detector. The data we use 
are the AURIGA filtered data obtained by means of the 
Wiener filter matched to a S-\ike signal, without setting 
any threshold on their amplitudes. If a signal enters the 
detector at time to, its output can be written as 



y(t) = hf(t-t ) + r)(t) , 



(2) 



where f(t — to) is the normalized signal template of the 
Wiener filter, h its amplitude (i.e. /(0) = 1), to its arrival 
time and T)(t) is a stochastic process with zero mean and 
correlation 



( V (t), V (t'))=a 2 h f(t-t'). 



(3) 



Here a\ is the variance of the filter output in the absence 
of any signal and f(t) is a superposition of two expo- 
nentially damped oscillating functions p5[ | which can be 
approximately expressed as 

fit) w e — 1*1/*- cos(w i) cos(u> s i) (u t w > 1) , (4) 

where t w is the Wiener filter decaying time (i.e. the in- 
verse of the detector bandwidth) , loq a center carrier fre- 
quency and ujb is an amplitude modulation frequency. 
Typical values for the AURIGA detector are t w ~ 1 s, 
ujo ^ 920 Hz and uj b - 20 Hz. 

Let us define now the random variable X, which rep- 
resents an averaged measure of the energy released by a 
g.w. signal impinging on the detector at time to 
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2W 
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dt\y(t~t )\ 2 



(5) 
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where t is the center of the time window. If an associ- 
ation between GRB and gravitational waves exists, the 
filtered output of the gravitational waves detector, in pe- 
riods just prior the GRB ("on-soured population) will 
differ (statistically) from the output at other times ( "off- 
source" population). 



A statistically significant difference between on- and 
off-source populations clearly supports a GWB-GRB as- 
sociation. For each of the N on GRB trigger we compute 
X(t^), k — 1, . . . , N OD which forms the Xon set of on- 
source events, and construct a complementary set Xoff 
with N a ff off-source events using windows before and af- 
ter the trigger. The sets Xon and Xoff are samples drawn 
from the populations whose distributions we denote p on 
and Pag. The off-source events are taken in periods not 
correlated with the trigger time, at a distance in time 
greater than 10 3 seconds from GRB the trigger, both be- 
fore and after it; this should be sufficient to have a fair 
sample of the off-source events as any GRB-GWB asso- 
ciation is reasonably excluded. 

For windows W greater than the Wiener filter charac- 
teristic time, the central limit theorem implies that p s 
is a normal distribution. 

Now suppose that the GWBs fall within the window 
W opened around the GRB trigger time and that the 
SNR of the gravitational signal associated with the GRB 
(averaged over the source population) is smaller than one, 
then p on is also a normal distribution with mean 



Hon = Moff + E 



1 



dt\hf(t-t 7 )\' 



liag+\^]E[h J ] (IU»U), (6) 

where E [ • ] is the average over the astrophysical source 
population of GRBs. 

The basic idea behind the statistical approach is a hy- 
pothesis testing where the null hypothesis Ho to test is 
the equivalence of the off-source and on-source distribu- 
tions: 



H Q : p ff(X) = p on (X). 



(7) 



The rejection of Ho clearly supports a GWB-GRB associ- 
ation. Since p on and p g are normal and could differ only 
in their mean values, we can test Ho by the Student's 
t-test |||. 

The t statistic is defined from Xon and Xoff by 



t = 



Hon - Hoff / N on N, 



v /r 



N OD + N off 



(N on - l)a 2 on + {N off - 1)<7 



oir 



N on + N off -2 



(8) 
(9) 



where pL n and £t ff {&1n and cr^g) are the sample means 
(variances) of Xon and Xoff, respectively. 

The expected value of t averaged on the source popu- 
lation and on the filtered output of the detector is 



ft=E [t] = f^gJj^, (10) 



where a — E[S]. 
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Let us consider the upper limit on E[h 2 ) in the assump- 
tion H.Q is true: in this case the most probable value of 
E[/i 2 ] is zero. From Eq. [l(] we get 



(tvA E[^] < 
\8W 



N off 



N nn N, 



oil 



f^t^max [N m = N off = Ny) 

flt,max 



(11) 



where fit.max is the upper limit. Assuming TLq true and 
N g > N on , we have 



V{h 2 } < h 2 max 



\t w ) JW 



(12) 



the value of fit,max can be deduced by the selected con- 
fidence level. In this way we are able to set an upper 
limit on the average amplitude of gravitational signals 
associated with GRBs 



1.96 



(N on 

x 

\100 



-1/2 



[5 x 10 



-1912 • 



(13) 



III. RESULTS 

The AURIGA data used in the two analysis are relative 
to the years 1997 and 1998, and the number of GRBs 
which fall into the AURIGA data taking periods is 120. 



A. Coincidence search 

Within the window of 5 seconds we have found 2 events 
in coincidence. This experimental result has to be com- 
pared with the number of coincidences due to chance. 

The shifts method consists of 10 4 time shifts of the co- 
incidence window; for each one we compute the number 
of coincidences; if the candidate GWB arrival times and 
the GRB triggers can both be modeled as Poisson ran- 
dom points p7[ , the number of accidental coincidences 
is fitted to a Poisson curve; from the fit we get the ex- 
pected number of coincidences due to chance. The re- 
sults are summarized in Fig. From the fit we ob- 
tain a value of mean expected accidental coincidences of 
n a = 2.57 ±0.04. 

Another approach to evaluate the number of acci- 
dental coincidences is given by Eq. [IJ that holds in 
case of Poisson random points ^?J. The total number 
of AURIGA events and GRB triggers are respectively 
N A = 26816, N 7 = 120, assuming AT = 2W = 10 s and 
T = 1.32 x 10 7 s we get n a = 2.4 ± 0.2. The error on n a 
can be easily estimated assuming the Poisson statistics 
for the fluctuations on Na and N~ i.e. 



respectively. The two estimates of the accidental number 
of coincidences are in good agreement and demonstrate 
that the two event rates are uncorrelated. We conclude 
that the 2 coincidence found are due to chance. 



B. Statistical search 

We have first tested the method using a Monte Carlo 
simulation by adding N on signals at fixed SNR over a 
gaussian noise generated by means of a noise model with 
the same parameters of the AURIGA detector. The sim- 
ulated detector output is then fed to the same data filter- 
ing procedures used for the AURIGA experimental data. 
We have generated 120 <5-like signals with SNR 3, 2 and 
1 and 1000 with SNR 1 storing their true arrival times 
to. The signals have been superimposed to stationary 
gaussian noise and we have then formed the on and off 
populations and calculated the t value. The probability 
P{t) that the t value obtained is due to chance is reported 
in Table § For signal with SNR = 3 and SNR = 2 this 
probability is very small and the value of t is statisti- 
cally significant, showing a GRB-GWB association. On 
the other hand, for signals with SNR — 1, we have to 
increase the number of GRBs to 1000 to get a statistical 
significance. 

The XoS and Xon sets relative to the AURIGA-BATSE 
data between 1997 and 1998 are given in the histograms 
of Fig. U and Fig. |3[ where the non gaussian tails are due 
to the non stationarity of the AURIGA noise. The value 
of the Student's t-test obtained is t = 0.58, which corre- 
sponds to a probability of 0.28 that it is due to chance. 

As we conclude that there is no evidence of an asso- 
ciation of GRB-GWB, we can put a constraint on gra- 
vitational radiation emitted from GRBs averaged over 
the source population, /irms- The value of fit.max can 
be found setting a confidence level C.L. and solving the 
equation: 



f d (t) dt = C.L. 



(14) 



where fd(t) is the distribution function for the Student's 
t with d degrees of freedom. Notice that for high d, fd(t) 
tends to a normal curve with zero mean and unitary vari- 
ance. If we choose C.L. =95% and set d = N on + N a g — 2, 
we get fit, max = 1-65. Using Eq. [TJ], we set the following 
upper limit on the averaged g.w. amplitude associated 
with GRBs: 



'RMS 



< 1.5 x 10" 



(15) 



IV. DISCUSSION 

The two search methods reported in this paper show no 
evidence of a correlation between 7-ray bursts and gravi- 
tational waves: First, the two coincidences found with 
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the coincidence search have no statistical significance 
and can be reconduced to chance. Next, the statistical 
method doesn't lead to a statistically significative value 
of the Student's t. However we were able to put an up- 
per limit on gravitational signals associated to the GRBs 
averaged over the source population, /irms < 1-5 x 1CP 18 . 

The existence of burst-like excitations (usually referred 
as non-gaussian noise) in the AURIGA data is well 
known, and here we deal with this noise selecting the 
periods of time when the detector was operating in a sat- 
isfactory way. The vetoing procedure of data dominated 
by non gaussian noise has been set by the AURIGA data 
analysis ||8[ and the resulting duty cycle is about 40 % 
of the total operating time. Care was also taken to cope 
with non stationarities of the AURIGA noise that give 
rise to the non gaussian tails in Fig. || and Fig. ^ ex- 
cluding such tails from the fits we used to estimate the 
statistical variable t. However, it is important to no- 
tice that even with the above limitations the analysis 
has reached a level of sensitivity which is astrophysically 
of some interest. In order to get better upper limits we 
have also explored the possibility of using in our analysis 
the incoming direction of GRBs but we have found no 
significant improvements (see Appendix |b|) . 

To increase the confidence in our results we have ap- 
plied the non-parametric Mann- Whitney it-test [E9| to 
the sample sets of the statistical search, obtaining a sec- 
ond confirmation of the null hypothesis. Ranking the 
elements of the union set Xoff © Xon m increasing order, 
we get a statistical parameter z = 1.59, smaller than the 
critic value z = 1.95 imposed by the fixed C.L. of 95 %. 

The upper limit we have set can be improved in the 
future simply by the increasing of common data taking 
periods of AURIGA (N on increases) and the new exper- 
iment HETE-II (High Energy Transient Explorer (30|) 
which is going to substitute by now the wasted BATSE 
satellite in the GRB search. Another possibility to reach 
more astrophysically interesting sensitivities is the up- 

§de of the AURIGA detector which is now in progress 
The predicted sensitivity and bandwidth of the 
AURIGA detector equipped with the new read out sys- 
tem would be respectively h m i n » 8 x 1CP 20 and i" 1 » 10 
Hz. This sensitivity, together with an enhancement of 
noise stationarity and duty cycle of the detector, would 
correspond to the lowering of the upper limit /irms of 
about 2 order of magnitude in one year of correlation 
analysis. 



P({k a ,t a }) = e- nt «/ T 



(nt a /T) k * 

k a l 



(Al) 



For to not overlapping windows, it can be demonstrated 
that in the limit of n — > oo, T — > oo and n/T constant, 
the probability of {fci points in t\}, . . . , {k m points in 
t m } is O 



P({h,h},...,{k m ,t m }) =]Jt 



_ nU/T {nU/Tf 



= ]JP({k i ,t i }) (A2) 



showing that the events {ki points in t{} and {kj points 
in tj} are independent for every i and j. Substituting 
m = N lt U = 2W = AT and n = N A we get 



P({k 1: AT},...,{k N ^AT}) 

k N y 



-N. 



N A AT/T ( N A AT ^ J_ 



n 

v ' i=l 



(A3) 



We can now evaluate the probability P(k) to have k 
coincidences, that can be obtained taking into account, 
at fixed k, all the possible sets [ki] = {k\, • • • , fc/v } with 



the constrain ki — k; we get 

P(k) = J2P AT}, AT}) 

[ki] 

k N-, 



-N- 



N A AT/T ( N A AT \ 



\ T 



\ki\ i=i 



-N-,N A AT/T I N A AT\ ' iV* 



T 



k\ 



(A4) 



Rearranging the factors in Eq. A4 we get a Poisson dis- 
tribution as in Eq. [I]. The last equality can be easily 
demonstrated using the multinomial expansion relation 
1321 



N 



5> = k ^W±j 



(A5) 



\i=l / [ki] «=1 

and substituting n = 1, i = 1 • • • N. 



APPENDIX A: 

To derive Eq. ^ let us consider n points random 
distributed in a time interval [0, T]. The probability 
P ({k a , t a }) that k a points lies in the time window t a — 
ti — ti is given by the binomial distribution |2j] ] with 
probability p = t a /T that a single point lie in t a . If 
n>l and t a T, using the Poisson theorem we get 



APPENDIX B: 

As the incoming direction of GRBs is known, one may 
wonder if a selection of the GRBs based on a cutoff on 
the AURIGA antenna pattern, averaged on the unknown 
polarizations of the GWB, could increase the sensitivity 
of our analysis. 

The antenna pattern of a resonant bar detector is given 

by 
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F{6) = 1- (k-z^) 2 = sin 2 (6>) 



(Bl) 



where fc and z are unitary vectors parallel to the GRB 
direction and the antenna bar axis and 6 is the angle 
between k and z. Sources that fall outside the two cones, 
which arc defined by the equation k • z > cos(£), have a 
figure pattern F(9) > sin 2 (£) and therefore the average 
energy associated with these g.w. sources is 



75 
G4 



sin(0 + 7* sin(3£) 
6 



50 



sin(5£) 



(B2) 



where the solid angle fij is defined by F{&) > F{£). 
Moreover, as GRBs are isotropically distributed over the 
sky, the cutoff on the figure pattern decreases the num- 
ber of available GRB: 7V|„ = iV „sin(£). Therefore the 
net effect of this cutoff on the expected value of t in Eq. 
Col is 



75 



1 



G4 

x sin(£) 1/a 



sin(£) + - sin(3£) 



G 



1 

50 



sin(5£) 



(B3) 



The function nf/^t is a continuously increasing function 
in the range £ G [0, 7r/2], and n\/nt = 1 at £ = tt/2 (i.e. 
the whole solid angle). We must conclude that a cutoff 
on the GRB direction does not enhance the sensitivity of 
the statistical search. 
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TABLE I. Results of a Monte Carlo simulation of the sta- 
tistical search. The column Pit) is the probability that the 
estimated t is due to chance. 
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N on 


t 


P(t) 
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120 
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120 


3.4 


4 x 10~ 4 
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0.4 
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3.6 


KT 4 
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